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I. INTRODUCTION 


Constituent quark models, on one side, and the parton picture, on the other side, repre¬ 
sent two complementary descriptions of hadron structure |l|. The rest frame (or laboratory) 
description, based on constituent quarks, and the infinite momentum frame description 
(Poo), based on partons, appeared as formalisms capable of interpreting particle physics 
phenomena at different resolution scales. The algebraic approach, we are going to explore, 
seemed useful in establishing links between them [5jj. It was intensively investigated in the 
late sixties, but a quantitative scheme for interacting scenarios was never attained. 

The birth of the Quantum Chromodynamics (QCD) set the general framework to under¬ 
stand deep inelastic phenomena beyond the parton model. The perturbative approach to 
QCD is able to connect observables at different resolution scales, but the realization of the 
complete project, i.e., to fully understand the consequences of the dynamics of quarks and 
gluons, requires the input of unknown non-perturbative matrix elements to provide absolute 
values for the observables at any scale. 

In the recent past the work by Gluck, Reya and coworkers || has shown that, starting 
from a parton parametrization at a low resolution scale /i q, the experimental deep inelastic 
structure functions at high momentum transfer can be reproduced and even predicted [£ij. 
/i ( 2 , is evaluated by evolving back the second moment of the valence distribution to the point 
where it becomes dominant. The procedure, suggested by Parisi and Petronzio ||, assumes 
that there exist a scale, /i 2 , where the short range (perturbative) part of the interaction is 
negligible, therefore the glue and sea are suppressed, and the long range (confining) part 
of the interaction produces a proton composed of three (valence) quarks only. Jaffe and 
Ross [|] proposed thereafter, to adscribe the quark model calculations of matrix elements 
to that, hadronic scale, /i ( 2 . For larger Q 2 their Wilson coefficients will give the evolution 
as dictated by perturbative QCD. In this way quark models, summarizing a great deal of 
hadronic properties, may substitute for low-energy parametrizations. 

Following such a path, a partonic description can be generated from gluon radiation even 
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off a purely valence quark system, which can be used to generate the non perturbative input 
occurring in the Operator Product Expansion (OPE) analysis of lepton-hadron scattering 
in QCD [j7j. A systematic analysis based on non-relativistic potential models, shows that 
the approach can be consistently developed at Next-to-Leading Order both for polarized 
and unpolarized structure functions [Q, including non-perturbative contributions from the 
nucleon cloud || or from the partonic substructure of the constituent quarks [dOi]. The 
formulation of this naive approach has been motivated by simplicity and by the success of 
the constituent models in reproducing many properties of the hadronic spectrum in the 1 - 
2 GeV region. 

It has been quite evident, since the original formulation of many of these models, that rel¬ 
ativistic effects to the nucleon wave function, as well as, covariance requirements, are needed 
even for a phenomenological description of the structure of hadrons. The present paper is 
devoted to the study of a possible generalization of our approach in that direction, and we 
will show that relativistic covariance can be incorporated within the same framework in a 
rather transparent and elegant manner. To this aim we develop a constituent quark model 
in the light-front realization of the Hamiltonian dynamics and apply to it the formalism for 
the calculation of both polarized and unpolarized parton distributions. 

In the next section we proceed to reformulate the procedure using the Hamiltonian 
formalism in the light front dynamics. A light-front constituent quark model is defined and 
its momentum distributions are related to the quark distribution function at the hadronic 
scale in the conventional way. The third section is dedicated to analyze, under different 
hadronic conditions, the perturbative evolution required to reach the scale of the data. Two 
scenarios, characterized by the presence or absence of soft gluons at the hadronic scale, are 
studied. We end by summarizing our results and discussing those points, which at the light 
of our presentation, require further development. 
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II. QUARKS AND PARTONS 


A schematic phenomenological distinction between constituent quarks and partons is 
often made by differentiating the infinite momentum frame, where partons show up, from 
the laboratory rest frame, where the constituent quarks appear to be the relevant degrees 
of freedom responsible for the correct symmetries and quantum numbers. Of course such a 
schematic point of view is only approximate and can also be at the origin of some misun¬ 
derstanding. Actually the two descriptions appear well defined in the light-front description 
of deep inelastic scattering, where the parton model is recovered, in the Bjorken limit, due 
to the dominance of short light-cone distances ( z 2 ~ 0) in the relevant Feynmann diagrams. 
As a consequence the partonic description can be developed in the rest frame of the hadron 
by using light-cone formalism. In particular the i-th parton distribution can be related to 
the light-cone momentum density^. 


<?i IU) 0,ho) = 


(■ 1-xf 


d 3 k 4 (i) (k 2 )<5 


x 


k + 


1 — x Mi 


N , 


( 1 ) 


where k + /P + = k + /Mjy = (^/k 2 + mj + k z )/M^r is the light-cone momentum fraction of the 
struck parton in the rest frame, Mn and m,; are the nucleon and parton mass respectively 
and nj (k 2 ), r4(k 2 ) represent the light-cone density momentum density of the i-th parton 
whose spin is aligned (j) or anti-aligned (j) to the total spin of the parent nucleon. If one 
assumes that at the scale /j ( 2 only the u and d constituent quarks are resolved, the momentum 
densities can be written 


d(l)(k 2 ) = (iV,J« = +l/2|X: 

2=1 

The light-cone distributions (|2|) have been evaluated in the past making use of non-relativistic 
constituent quark models, while in the present investigation we want to improve their descrip¬ 
tion including relativistic effects as introduced by a light-front formulation of a three-body 


t? 1 + (~)a z 

— L 1 S(k-k t )\N,J z = +l/2) . (2) 


X A formal derivation of Eq. (]lj) can be found in ref. |8j. Eq. (!]) includes the support correction 
as discussed in the same reference. 
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interacting system. As a consequence we will remain within a constituent picture where 
the partons in the rest frame are identified with three (constituent) quarks at the hadronic 
scale, and covariance requirement is fulfilled. 

Since the hadronic scale //q turns out to be very low (//§ ~ (0.1 — 0.2) GeV 2 ), close to 
the constituent quark massQ we assume that the constituent picture at this scale has to be 
identified with a Constituent Quark Model, whose parameters are fixed to reproduce the 
basic features of the nucleon spectrum in the energy region 1 — 2 GeV. 


A. The Light-Front Constituent Quark Model 


Before discussing the details of the quark wave functions, it is worth to remark the main 
features of the theoretical framework that we have employed to achieve a covariant descrip¬ 
tion of the nucleon’s constituent degrees of freedom. An exhaustive review of covariant 


hamiltonian dynamics can be found in [11 


On a quantum level, the covariance requirement is formulated by demanding the proba¬ 
bility assigned to a given event to be the same in all inertial frames. According to Wigner 
theorem El. this is achieved if (and only if) the correspondence between the states in 
different frames is realized by means of an unitary representation U (A, a) of the Poincare 
group. 

The connected component of the Poincare group is a ten parameter Lie Group that 
includes four space time translations, three rotations and three Lorentz boosts. Let us 
denote with P M the representations of the generators of space-time translations, with K 2 
those of Lorentz boosts, and with J J those of SU( 2). Any linear combination of generators is 
still a generator. To achieve a covariant formulation of the Poincare Algebra it is convenient 
to introduce a new set of tensor generators of the Lorentz group: 


2 The actual value of the scale /j,q ranges from 0.094 GeV 2 , if only valence quarks are considered, 
to 0.37 GeV 2 @, when non-perturbative q — q pairs and gluons are included. 
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J 0j = K j 

(3) 

jjk = pkij i 

(4) 

jap _ _jP& 

(5) 


The tensor character of these operators is expressed by the following relationships: 


f/ t (A)P M f/(A) = A P P U (6) 

U\A)J pv U(A) = A P A v a J pa . (7) 

Starting from the equations © and © it is possible to work out the following Algebra : 

[P^P V ] = 0, (8) 

[J pp , = i^pP _ (9) 


[J p \ Jp°] = i{gt*r p - g pp J ua - g ua J pp + g up J pa ). 


( 10 ) 


Hence, in order to build a covariant quantum mechanical model, one has to construct a 
model Hilbert space and find a set of operators such that these commutation relations are 
satisfied. 

In general, such a representation is not unique and depends on the interaction. However, 
it is possible to single out a subset of generators which do not contain the interaction. Such 
operators group together to form the kinetic subgroup. Once the kinetic subgroup is deter¬ 
mined, the whole set of generators is also fixed including the Hamiltonian generators (i.e. 
those operators which contain the interaction). Dirac proposed [H| to fix first the kinetic 
subgroup identified with the stability group of a three-dimensional surface in Minkowski’s 
space, and then extract the Hamiltonian generators by group representation properties. 

The choice of the surface that determines the kinetic subgroup fixes the form of the 
dynamics. In the following we will use the light-front form of the dynamics whose invariant 
surface is the light-front (null) plane x° + x 3 = 0. This form of the dynamics has many 
important properties, for example, it leads to the largest possible kinetic subgroup. 
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In the light-front form of the dynamics, it is convenient to consider the following set 
of generators: P 1 , P 2 , P + = P° + P 3 , J 3 , K 3 (kinetic subgroup) and P~ = P° — P 3 , 

= Ki + z x Jj_ ( set of Hamiltonians). In particular, the role of the Hamiltonian H 
of non relativistic Quantum Mechanics is played, in the front form, by the generator P~ of 
translations along the direction orthogonal to the null plane. This operator has a positive 
eigenvalued spectrum and this prevents the existence of negative ” energy” states, allowing 
for the number of particles to be fixed. For this reason the light-front form of the dynamics 
is particularly well suited to describe a low energy model where only three valence quarks 
are considered. 

The construction of a representation of the Poincare space such that the kinetic subgroup 
is the stability group on the light-front, is clearly much easier by using a parametrization of 
the Minkowski space which is as close as possible to the symmetry of the problem, i.e. the 
light-front coordinates: x + = x° + x 3 , x~ = x° — x 3 , a; 1 and x 2 . In this way a vector is given 
by (x + , x~, x 1 , x 2 ) = (x + , x~, xj_), an the scalar product bjx-y — ~\x + y~ — \x~y + +^.±- yj_. 
The null plane is the plane with x + = 0 and the ” (three)-vector” (i“,xi) on that plane is a 
true three-vector under light-front boosts (i.e. the boost belonging to the kinetic subgroup). 

In order to model the Hilbert space one needs to combine the generators of the Poincare 
group to get a set of self commuting operators which have a physical interpretation. By 
combining P M , J pv one can construct two invariants: 

i) the mass operator 

M 2 = -P 2 = P+P” - P 2 ± , 

ii) and the square (HQ W p ) of the Pauli-Lubansky vector: 

W = - P u J pa e 

vy /x 2 u • 

In relativistic systems the definition of the spin operator is related to the Pauli-Lubansky 
vector: one defines, in fact, the square of the spin by 

W 2 = —M 2 j 2 . ( 11 ) 
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It can be proved that there is an infinite set of operator-valued vectors satisfying the relation 
©• Therefore, the spin operator is defined as a linear combination of and M h 


Ja = K(P)W,M 


-1 


( 12 ) 


The coefficients u^{P) are space-like operator-valued vectors of an orthonormal basis |T4 . 
Different choices of this basis leads to the specific definition of spin in different forms of the 


dynamics. Starting from fll2|) one can prove that the spin operator transforms as 


UH.\)]U(\) = R W (.\,P)J, 


(13) 


where Rw(A, P) a p = U a (AP) ■ A up(P) is called Wigner rotation. One can demonstrate that 
defining u 0 (P ) = P/M, the 16 components quantity u"(P) forms a -SO(1, 3) representation 
of a Lorentz transformation. The spin operator can be written 

1 1 n 

M = 5 eiw), «J(P) K(P) J f ° = yo« o,«J(P) <(P) £ If . (14) 

where the last equality holds when the system is made up of n independent constituents. 
In the instant form of the dynamics (in which the invariant surface is chosen to be the 
equal-time plane) the Wigner transformations associated to a rotation is the rotation itself 
and the spin is called canonical. One can demonstrate that, in this form of the dynamics, 
the ordinary composition rules of the angular momentum hold, providing one rotates the 
constituent spins by means of appropriate Wigner rotations. In the case of light-front form 
of the dynamics, the Wigner rotations associated to a light-front boosts are the identity. The 
ordinary composition laws for the total angular momentum hold, provided we transform the 
light-front spin into the canonical spin. The transformation that relates the light front spin 
to the canonical spin is called a Melosh rotation, 


Rm(P)c0 ■■= 


Its spin-1 representation is: 


[R m {P,M)] 


M + P + — io ■ (z x Pj_) 


(15) 


(m + p+y + Pi 


(16) 








1. Non-interacting and interacting systems on the light-front 


The single spin-1/2 particle states can be described by wave functions 0(p, A), where 
p = ( p + , p_|_) is the momentum three-vector on the light-front and A the longitudinal spin 
component. One has: 

m 2 = e7#wj> + ,pji 2 . 

A =-j J P 

d 3 \p] = dp + d 2 p ± 6{jp + ); (17) 


where the measure d 3 [p]/p + and the component A are invariant under the kinematic trans¬ 
formations, and the factor l/VE can be included in the definition of the wave function 
provided we change its light-front boost transformation properties accordingly. 

For a many-body system of n non-interacting particles the wave function can be written 
$(pi,Ai;...;p n , A n ), and Eq. flT7|) becomes 


l$ll 2 = E 

Al 


ST) 


|4(pi,Ai;...;p„, A„)|-’ . 


(18) 


Quite often it is more useful to represent <f> as functions of the total momentum on the 
light-front P = Ya=\ Pi: th e momentum fraction Xi = pf / P + , and the transverse relative 
momenta kjj_ = u± ■ Pi = pjj_ — XjPj^. Eq. (|18|) can be then transformed into 


l$ll 2 = E 

Al ,...,A n 


d 3 [P} 
P+ 




vi=l X r 


\i = 1 


Ki =1 


x <5(-P, k 1 ± ,xi,Ai; \ r 


(19) 


One can introduce a further set of coordinates adding, on kinematical basis only, a third 
component k {3 to the transverse momenta k.jj_ boosting the lab momenta pp. 


ki3 = u 3 (P) -Pi= - 
UJi = -Uo(P) -Pi=\ 


M 0 x t + 

M 0 Xi 


= M 0 Xi - UJi , 


M 0 Xi 


+ kf x 


M 0 Xi 


( 20 ) 


where the total mass operator of the system 
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( 21 ) 


Mo = 5] y p? + rnf = 


\ 


E 


k?i + mf 


Xi 


- E ^ — E v k ’ ? + 


m? . 


The invariant measure to be used in the definition of the scalar product transforms in the 
following way: 




U=] ^ 


i_L 




<i=l 


li k ; 


V.i=l u i 


The total spin, in the new coordinate system, is 


j = E [ l W , Xk,+ 

1=1 


Si 


with 


( 22 ) 


(23) 


s, = R M (ki,m,) ji (24) 

and the representation of the Melosh rotation reads 

n i r d n u + x i M o ~ ' (* x k i± ) 

[^(ki.rrii)] =- -- • (25) 

V C m t + XiM 0 ) 2 + k- L 

As a first consequence of this section we can conclude that the dynamics of n non¬ 
interacting spin 1/2 particles can be described, on the light-front, by a wave function 
-0(P, k 1; fj,i, k n , n n ) whose norm is given by 


= (?/#> = E 


d 3 P 


. p+ , 

\i= 1 


11 1- ! s IZ k .) W' , 


(26) 




where Hi are the eigenvalues of the spin projections s\ of eq. (p4|). We stress that the Jacobian 
and the invariant measure of the transformation (kjj_, —> k,, has been absorbed in the 
definition of the wave functions -0 as can be seen from Eq. (j26j) (cfr. the discussion after 
Eq. (0). 

The extension to interacting systems requires a dynamical representation of the Poincare 
group. One way to achieve this is to add an interaction V to the free mass operator Mq 


M — M 0 + V , 


(27) 
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where V is the interaction. All the other definitions remain unaffected, including the def¬ 
inition of the third component of the relative momenta ki 3, because the boosts we use are 
interaction independent. All required commutation relations are satisfied if the mass opera¬ 
tor commutes with the total spin j and with the kinematic generators. In the representation 
in which states are represented by functions of P, k, and /q these conditions are realized if: 

i) V is independent on the total momentum P; 

ii) P is invariant under ordinary rotations. 

Summarizing: in the light-front formulation of the quark dynamics, the intrinsic momenta 
of the constituent quarks (ki) can be obtained from the corresponding momenta (pt) in a 
generic reference frame through a light-front boost (ki = u(P ) • pi, P = X)f=i Pi) such that 
the Wigner rotations reduce to the identity. The spin and spatial degrees of freedom are 
described by the wave function: 

v b = ~^=5(P -p)x( ki,hi,---,k 3 ,/i 3 ), ( 28 ) 

where /q refer to the eigenvalue of the light-front spin, so that the spin part of the wave func¬ 
tion is transformed by the tensor product of three independent Melosh rotations: P,^ (k,, rrii) 
@ , namely V) = n|,;=i ^(k,, rrq). 

The internal wave function is an eigenstate of the baryon mass operator M = Mq + V, 
with M 0 = X)i =1 = X)i=i \J kf + rnf , where the interaction term V must be independent 

on the total momentum P and invariant under rotations. 

The nucleon state is then characterized by isospin (and its third component), parity, 
light-front (non-interacting) angular momentum operators J and projection J n , where the 
unitary vector n = (0, 0,1) defines the spin quantization axis. 

B. Valence Quark Hamiltonian 

I11 the present work we will discuss results of a confining mass equation (|27|) of the 
following kind 
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(Mo + V) iMO = (E H + m i - J + K ^) = MV>o,o(0 , (29) 

where k* = 0, £ = \J p 2 + A 2 is the radius of the hyper-sphere in six dimension and p 
and A are the intrinsic Jacobi coordinates p = (r 1 — r 2 )/v / 2, A = (ri + r 2 — 2r 3 )/ v / 6. The 
choice of the mass operator (^) has a twofold motivation: 

i) Simplicity. This is a first attempt to develop a covariant approach to DIS based on quark 
models, therefore we use a mass operator whose symmetry properties facilitate the numerical 
solutions. In this respect the choice of a hypercentral potential has big advantages. 

ii) The nucleon spectrum is reproduced. This has been demonstrated both within non- 
relativistic |H| and relativistic 13 approaches []. In particular the well known problem of 
the mass of the Roper is solved in the present case by the use of r/£ potential, as discussed 
by Ferraris et al. |L5|. Of course our aim is much less ambitious than reproducing all the 
complexity of the baryon spectrum, we restrict the calculation to the nucleon wave function 
and we do not consider spin-dependent terms assuming pure SU (6)-symmet.ric stated 

The intrinsic nucleon state is antisymmetric in the color degrees of freedom and sym¬ 
metric with respect the orbital, spin and flavor coordinates. In particular, disregarding the 
color part, one can write 


3 In ref. (l(| a simplified version of the same mass operator has been introduced to investigate 
the baryonic mass spectrum in the instant form dynamics. Namely the eigenstates of the operator 
M 2 = + nr 2 ) — a/£ + b £ are discussed, and the potential — a/£ 4- b £ is the simplified form 

of the more complex expression V 2 + {Mo, V} + XJ/j y k 2 + m 2 k 2 + m 2 obtained by squaring 
the correct mass operator M = Mo + V. We will discuss solutions of the mass operator Eq. ( |29|) 
directly, without further simplification. 

4 Adding a perturbative hyperfine interaction as discussed in ref. [15] in the context of non- 
relativistic Hamiltonians, would be rather simple also within the relativized scheme of Eq. (|29|). 
At the present stage we neglect such a complication, but a more complete study of the nucleon 
spectrum should include S'I7(6)-breaking terms within the light-front approach. 
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( 30 ) 


\at T T i / o \ / / -I > ( 0 , 0 ) / n \ Xms&ms + Xma^ma 

\N, J, Jn = +1/2) = V+o(0 y [ 0 , 0 , 0 ] (W -- 

where ^ 7 ,„(£) is the hyper-radial wave function solution of Eq. ([29]), y y (^) the hyper- 
spherical harmonics defined in the hyper-sphere of unitary radius, and 0 and x the flavor 
and spin wave function of mixed SU( 2) symmetry. Let us note that, in order to preserve 
relativistic covariance, the spin wave functions 


Xms = -j= [2 TU - (U + IT) t] ; 


xma = (U - IT) T , 


(31) 


have to be formulated by means of the appropriate Melosh transformation of the i-th quark 
spin wave function: 


Ti= Ru^mi) 


li 


f A/r \ 

rrii + Xi M 0 

1 

w 

\J ( rrii + XiM 0 ) 2 + k? ± 

y kiR J 


(32) 


U= RM(K,rrii ) 


M 

1 

( _ k . T \ 

oJ 

\j (rrii + XiMo ) 2 + kf x 

K mi + Xi M 0 


(33) 


where x t = pf / P + = (k iz + uy)/M 0 and k L/R = k x ± i k y . 


1. Numerical Solutions 


We have solved the mass equation (|29|) numerically by expanding the hyper-radial wave 
functions ^ 7 „(£) on a truncated set of hyper-harmonic oscillator basis states. The matrix 
elements of the mass operator (P9[ ) have been calculated in momentum space for the ki¬ 
netic energy term, and in configuration space for the interaction one. Making use of the 
Raylcigh-Ritz variational principle the HO constant has been determined and convergence 
has been reached considering a basis as large as 17 HO components [[Hj]. The parameters of 
the interaction, have been determined phenomenologically in order to reproduce the basic 
features of the (non strange) baryonic spectrum up to ~ 1500 MeV, namely the position of 
the Roper resonance and the average value of the 1“ states. We obtain: r = 3.3 and /q = 1.8 
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fm 1 [jTj to be compared with the corresponding non-relativistic fit r = 4.59 and Ki = 1.61 
fm 2 |L5|. The constituent quark masses have been chosen m u = m c i = m q = Mjy/3. 


The result of this reparametrization is that a huge amount of high momentum compo¬ 
nents is generated by solving the mass equation (cfr.Fig. la), and they play an important 
role in the evaluation of transition and elastic form factors within light-front constituent 


quark models as discussed by Caldarelli et al. [|18|] in connection with the solutions of the 
Isgur-Capstick model Hamiltonian. From the point of view of DIS the enhancement of high 
momentum components in the density distribution of the valence quarks is crucial to re¬ 
produce the behavior of the unpolarized structure functions for large value of the Bjorken 
variable x as it will be discussed in the next section where also polarized responses are 
investigated at the hadronic scale. 


C. Partons at the Hadronic Scale 


The calculation of the partonic content at the hadronic scale as given by Eq. (0), is rather 
involved mainly because of the spin dynamics accounted for by the Melosh rotations (|32|) 
and (|33|), but finally the polarized distributions at the hadronic scale can be written in a 
rather elegant way: 

a Q (m + Jk'i + t + m 2 + k z ) +b Q t 

,t) -E- J - D(b,t) (34) 

(m + \jk 2 + t + m 2 \ +t 


o) = 


7r 


M 


N 


9 (1 — x ) 2 Jo 


dtn(k 


where Q = £ = kj_, 


D(k z ,t) 


\Jt + k 2 z + m 2 

yft -\~k 2 z + m 2 + k z 


(35) 


x (t + m 2 ) (1 — x) 

1 — x Mjj x 

and n[k 2 , k^_) is the total momentum density distribution of the valence quarks in the nucleon 
calculated making use of the eigenfunction i/j 00 of Eq. (|29|): 


k z (x, t) = 


M; 


N 
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n( k 2 ) — 3 /|* 0 | J «(k- k 3 ) 3('y k, l d 3 k!d 3 k 2 d 3 k 3 , 


( 36 ) 


with / d 3 kn(k 2 ) = 3. The coefficients ag and &g take the values a u r = b u i = 5; b d t = a rf i = 2; 
6 u t f • 


1. Effects of Melosh Rotations 

The effects of high momentum components on the unpolarized parton distributions 
J2 q =(u,d) o) + <? J (x, /id)) = u v{ x ,Ro) + dy(x, jil) at the hadronic scale are shown in 

Fig. lb. Their important role in reproducing the correct behavior of the structure functions 
for large values of the Bjorken variable x will be discussed in the next sections. The relevant 
effects of relativistic covariance are even more evident looking at the polarized distributions 
A uv(x,nl) = vl(x,n o) — Ad v (x,fil) = dJ(x,/j,l) — d^(x,n g)) where the spin dy¬ 

namics on the light-front plays a crucial role. The introduction of the Melosh rotations 
results in a substantial enhancement of the responses at large x and in an suppresion of the 
response for 0.1 < x < 0.5 as can be seen from Fig. 2. We show, in the same figure, also 
the predictions of a pure relativized solution obtained by solving numerically Eq. ( |29| ) and 
neglecting the Melosh rotation effects in (JTEI). Such a calculation retains the contribution 
due to the high momentum components, while the covariance requirement on the parton 
distribution is lost. One can argue on the quite relevant role of covariance on the spin 
observables. 


2. Non-perturbative Gluons 


A natural choice for the unpolarized gluon distribution within the present approach, has 
been discussed in refs. [|| and assumes a valence-like form (cfr. also ref. 0) 


ho) = [uv(x, nl) + d v (x, nl) 


(37) 
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This definition implies / G(x, /id) dx = 2 and therefore only 60% of the total momentum 
is carried by the valence quarks at the scale //q. In this case pi = 0.220 GeV 2 at NLO 
([«^(A i o)/( 47r )]NLO = 0.053). 

If the gluons were fully polarized at the scale pi one would have | AG(x, /x§)| = G(x, pi). 
However recent next-to-leading (NLO) studies of polarized Deep Inelastic Scattering (DIS) 
fl0| pO| show that the x-shape of AG is hardly constrained by the present data. As a 
consequence quite a few sets of parton distributions, with rather different :r-shapes of the 
gluon component, provide very good descriptions of the same data. In particular the sign 
of A G(x, Q 2 ) for x > 0.1 is not convincingly fixed: in the parametrization C of ref. |2l|] it is 
assumed to be negative , while it is taken positive in other cases (cfr., for example, ref. [[22||). 

In the recent past different authors have discussed the sign of the gluon contribution to 
the nucleon’s spin in the light-cone gauge ( A + = 0) and at the hadronic scale valid at the 
quark model level. In particular Ball, Forte and Ridolfi [Hj showed the compatibility of the 
experimental data with a large and positive gluon polarization at Q 2 = 1 GeV 2 (A G(Q 2 ) = 
Jq dxAG(x,Q 2 ) ~ 1.6 ± 0.9); Jaffe |26] showed that the color electric and magnetic fields, 
which are believed to be responsible for the spin splittings among light baryons, give rise to a 


significant negative contribution (AG(pl 


-0.7) at the hadronic scale; Mankiewicz, Pillcr 


and Saalfeld [^7[ performed a QCD sum rule calculation obtaining AG(pl) ~ 2.1 ± 1.0, 
while Barone, Calarco and Drago []28|, by using the non-relativistic quark model of Isgur- 
Karl and considering the renormalization due to self-interaction contributions, obtained 
A G(i4 = 0.25 GeV 2 ) « 0.24. 

We shall describe results in two scenarios characterized by different gluon distributions 
(AG) at the hadronic scale ||: 


i) Scenario A : A G(x, /ig) = 0. Only quark valence distributions are allowed at the 
hadronic scale. The momentum sum rule determines /i ( 2 = 0.094 GeV 2 at NLO 
([«s(ho)/( 47r )]NLO = 0.142). 

ii) Scenario B: A G(x,pl) = f G(x, p q). / is the fraction of polarized gluons and has to 
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be considered with the appropriate sign. For example the Jaffe’s suggestion, in our 
approximation scheme, would imply AG(x,nl) ~ — 0.35 G(x, /Xq), while / > 0 will 
imply positively polarized gluons Q 


III. RESULTS 


A. Partonic bremsstrahlung 


An important aspect of our approach is that we perform, in both the flavor non-singlet 
and singlet channels, a complete analysis at Next-to-Leading Order, an essential requirement 
needed to include the anomalous (75) contributions consistently. Complete NLO calculation 


of all polarized two-loop splitting functions |29|], in the MS renormalization and various 
factorization schemes |^4,^5,^9|), are now available. The evolution of the polarized distribu¬ 
tions is performed according to the solution of the renormalization group equation inverting 
the moments (f(Q 2 )) n = Jo dx f(x, Q 2 ) x n ~ l numerically. Since the starting point for the 
evolution (/ip), consistent with the renormalization scale of the quark model, is rather low, 
the form of the equations must guarantee complete symmetry for the evolution from /1 q 
to Q 2 / 1 q and back avoiding additional approximations associated with Taylor expan¬ 

sions and not with the genuine perturbative QCD expansion (a complete discussion on the 
evolution procedure can be found in ref. [El). 


We use the so called MS and AB factorization scheme evolutions on our model input to 
show that our results arc largely independent of which one we use. In principle the calculation 
should be independent of factorization scheme, however at present it is not straightforward 
to abscribe the non-perturbative input to any specific factorization scheme. The spurious 


5 In fact in ref. [|(J it has been shown that f AG(x, Hq) dx < 0. Such inequality does not imply 
AG(x,/ig) < 0 in the whole x-range, and our results for positively and negatvely polarized gluons 
can be seen as upper and lower limits for both cases. 
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factorization dependence arises because we associate the same non-perturbative input to the 
different factorization schemes as discussed in the next sections. 


B. Discussion 


1. Polarized Parton Distributions 


In Fig. 3 the results for the proton structure function gi(x, Q 2 ) are shown and compared 
with the experimental data. The non-relativistic approximation of the present calculation 
appears to reproduce rather poorly the experimental observations. Even the use of non- 
relativistic models which reproduce rather well the whole baryon spectrum does not alter 
this conclusion, as shown in ref. 0 where the predictions of quite a few non-relativistic 
models are discussed. 

The full covariant calculation leads to a theoretical predictions quite close to experimental 
data in the region 0.01 < x < 0.4 even within the simple assumption of a pure valence 
component at the hadronic scale (scenario A). The calculation is parameter-free and the 
only adjustable parameters of the Hamiltonian (r and Ki in Eq. (^)) have been fixed to 
reproduce the low-lying nucleon spectrum as discussed in section |11 B 1| . 

The effect of relativistic covariance in the quark wave function is mainly associated to 
the spin dynamics induced by the Melosli rotations (cfr. Eqs.(|32|), (fhil), (|54|)) and these 
transformations lead to a strong suppression of the structure function in the small-x region 
(x <0.5). 


In order to introduce the gluons non perturbatively we evolve the unpolarized distri¬ 
butions predicted by the scenario A, up to the scale of scenario B where 60% of the total 
momentum is carried by valence quarks. At that scale we add to the valence partons the 
non-perturbative gluons as defined in Eq. ( |i7|) of section |II C 2| . 

Moreover, at the new hadronic scale, the fraction of polarized gluons is chosen to be 


negative according to Jaffe’s result [26] or positive as in other parametrizations. Note that 


18 









by choosing a negative fraction we are maximizing the difference with respect to the radiative 
gluons, because those yield to a positive polarization. Fig. 6 leads us to confirm that the 
low- a; data on g\ do not constrain the gluon polarization strongly. If we vary the fraction of 
polarized gluons from 35% to 100% the quality of the agreement is not significantly changed 
in the region 0.01 < x < 0.4. For larger values of x the valence contribution plays a major 
role and the behavior of the structure functions might depend on the model wave functions. 

For both scenarios the boundary value /i ( 2 is rather low and one could question on higher 
order contributions to the evolution. In order to appreciate (and maximize) the relevance 
of NNLO (or higher order) effects, the polarized proton structure function are evolved at 
Q 2 = 3 GeV 2 within the scenario A and by using both LO and NLO approximations (cfr. 
Fig. 4). In the region of 0.01 < x < 0.5 they differ roughly of 20% which represents the order 
of magnitude of the largest higher leading effects. The same evolution obtained within the 
scenario B would lead to an even smaller difference (< 10%) between LO and NLO results. 

The comparison of the predicted neutron structure function with the data (Fig. 5) differs 
quite substantially according to the amount of polarized gluons at the hadronic scale. Within 
scenario A the values of gi(x,Q 2 ) remain quite small according to the fact that the mass 
operator chosen (|2U|) is SU(6) symmetric. The introduction of negative gluon polarization 
brings the predictions of the present relativistic quark model closer to the experimental 
observations at least in the x > 0.1 region; a larger negative fraction of gluon polarization is 
favored by the data as can be argued also from the results of positively polarized gluons as 
shown in the same figure. Any how, since the neutron structure function is largely dependent 
on the nucleon wave function model, our assumption on SU(6) symmetry is the largest source 
of uncertainty and a more detailed study of g™ should include more complex wave function 
models. Another source of uncertainty in our evaluation of the neutron stucture function is 
the lack, at the hadronic scale, of (polarized) sea contributions which can play even a major 
role for x < 0.3 |T0f . 


All these limitations, however, do not obscure the main aim of our study, namely the 
discussion on the effects due the relativistic requirement of the quark wave function. 
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Fig. 6 is devoted to the gluon distributions. Within scenario A the polarized gluon 
distributions remain positive at the experimental scale as a result of the NLO evolution. 
On the contrary a large amount of negative gluon polarization is predicted within scenario 
B if one assumes anti-aligned gluons at the hadronic scale. The distribution is largely 
dependent on the polarization fraction. The dotted line of Fig. 6a shows the consequence, 
in the deep inelastic regime, of the Jaffe’s calculation at low energy. The absence of non- 
perturbative sea polarization results in a huge (and probably unrealistic) amount of negative 
gluon polarization to reproduce the neutron data. Such a large amount of gluons is, however, 
inconsistent with the unpolarized gluon distributions as it is shown in Fig. 6b (dotted line 
again). 

Fig. 7 shows the results of the analysis on the factorization scheme dependence. We have 
chosen to show both the proton and the more appreciable case of the polarized neutron data. 
The small difference gives us confidence in the results. 

2. Unpolarized parton distributions 

Figs. 8 and 9 complete our study of unpolarized partons (gluons are shown in Fig. 6 as 
already mentioned). Particularly evident is the improvement, due to the high momentum 
components introduced by the relativized solution, in describing the proton F 2 distribution 
in the region 0.4 < x < 0.8. The non-relativistic description is largely suppressed for x > 0.4 
due to the lack of high momentum contributions. The solution of the mass operator (|29|) 
have the appropriate components to contribute in this x-region . The work is completed by 
the comparison of LO NLO evolution also in the case of unpolarized scattering (Fig. 8a) and 
the results on valence and sea distributions (Fig. 9). In this last case one could repeat the 
comments already formulated for the polarized distributions: the lack of non-perturbative 
sea at the hadronic scale largely limits the possibility of reproducing the parton distributions 
in the singlet sector. 
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IV. FURTHER DEVELOPMENTS AND CONCLUSIONS 


The present paper is devoted to a semi-quantitative analysis of the effects due to relativis¬ 
tic covariance requirement in the study of polarized structure functions within constituent 
quark models. The investigation clearly shows the relevance of relativistic effects due, in 
particular to the dynamics of the spin observables in special relativity. In fact two new in¬ 
gredients play a major role in our light-front calculation: i) the presence of high momentum 
components and ii) the spin dynamics induced by the Melosh transformations. 

i) The large amount of high momentum components in the nucleon wave function is gener¬ 
ated by the solution of the relativized mass equation (^) and are relevant in reproducing 
the large x behavior of polarized and unpolarized structure functions. The lack of such high 
momentum contributions have been repeatedly stressed in connection with non-relativistic 
calculations |8|,|30| without reaching a satisfying solution within non-relativistic quark mod¬ 
els. 

ii) The Melosh rotation dynamics introduce the basic new ingredient in the calculations and 
its effect is quite sizeable in suppressing the proton response in the region x ^ 0.4. 

The first natural extension of the present discussion is the study of transversity distri¬ 
butions as they will be investigated by Drcll-Yan processes. In fact relativity must play a 
quite sizeable role for that observable pT| [ki|| and the light-front dynamics is a privileged 
tool to study such effects. 

In addition to the advantages of our approach we want to stress also its limitations and 
approximations because a further study should remove them in order to provide a more 
quantitative study of the nucleon spin within the light-front dynamics. 

The most relevant limitation of the present investigation is probably the lack of (non- 
perturbative) sea contributions at the hadronic scale due to the meson cloud of the nucleon. 
The effects can be rather large and significant in particular for the spin response of the 
neutron. Work in this direction is in progress and we want to extend the investigation of 
non-relativistic approaches H| inserting the recent successful description of the meson cloud 
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34| within a coherent light-front formulation |35[ 


A second limitation is determined by the particular choice of the quark-quark interaction 
of Eq. (^). As mentioned in few points of the present paper the large x behavior of the 
nucleon response can be affected by the specific choice of the interquark interaction model. 
Work is in progress to investigate the numerical solutions of more sophisticated hyper-central 


potentials already discussed in the non-relativistic limit |15|. In addition we did not include 
hyperhne interaction terms in our analysis and we would like to extend the approach also 
to the consistent light-front treatment of such interaction. 

Finally the factorization scheme dependence is to be investigated in more detail in con¬ 
nection with radiative model studies of the kind discussed in the present paper. 
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FIGURES 




FIG. 1. Fig. la (left panel): the valence momentum distribution of Eq. (36) as function of 
|k|. The results of a full covariant light-front calculation (full curve) are compared with the 
non-relativistic approximation (dashed curve). The corresponding total valence distributions at 
the hadronic scale are shown on the right panel (Fig. lb). 




FIG. 2. Left panel (Fig. 2a): the polarized distribution Auy(x, //q) as function of x: the 


non-relativistic approximation (dot-dashed curve), and the relativized solution of Eq.(29) which 
neglects Melosh rotations (dotted curve) are compared with the results of a complete light-front 
calculation (full curve). On the right panel the distibution Ady (x, /xq) (same notation). 
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FIG. 3. The proton polarized structure function at Q 2 = 3 GeV 2 . The full curve represents 
the NLO (MS) results of a complete light-front calculation within a scenario where no gluons 
are present at the hadronic scale (scenario A); the corresponding non-relativistic calculation are 
shown by the dot-dashed line. Scenario B is summarized by the dotted line in the case of negative 
polarized gluon fraction (f AG = —0.7 as discussed in the text), and by the dashed line in the 
case of positive gluon polarization (f AG = +0.7). Data are from SMC [36], and SLAC(E142) [37] 
experiments. 
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FIG. 5. The NLO (MS) neutron polarized structure function at Q 2 = 3 GeV 2 within a complete 
light-front calculation. The full curve represents the results within a scenario where no gluons are 
present at the hadronic scale (scenario A). Scenario B is summarized by the dotted line in the 
case of negative polarized gluon fraction (f AG(x,/j,q ) = —0.7 as discussed in the text), and by 
the dashed line in the case of positive gluon polarization (/ AG(x,{Iq) = +0.7). Data are from 
SLAC(E154) [38] experiments. 
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FIG. 6. Left panel: Polarized gluon distributions at Q 2 = 3 GeV 2 obtained evolving at NLO 
(MS) the polarized partons in both scenarios. Scenario A (full line). Scenario B is summa¬ 
rized by the dotted line in the case of negative polarized gluon fraction (f AG(x,iiq) = —0.7 
as discussed in the text), and by the dot-dashed line in the case of positive gluon polarization 
(f AG(x,nl) = + 0 - 7 )- 

Right panel: unpolarized gluon distribution at Q 2 = 10 GeV 2 obtained evolving at NLO (DIS) 
unpolarized partons. Scenario A (full line); scenario B (dotted line). For comparison also the 
results at LO are shown in this case (dot-dashed). CETQ4 NLO (DIS) fit of ref.[39]: full triangles. 
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FIG. 7. Scheme dependence of the polarized proton and neutron structure functions at Q 2 = 3 
GeV 2 . MS dotted lines; AB: dashed lines. 



FIG. 8. Left panel: unpolarized F 2 structure function at Q 2 = 10 GeV 2 as function of x. LO 
evolution (dashed) NLO (DIS) full line. Right panel: the F 2 structure function of the proton at 
Q 2 = 10 GeV 2 as predicted by the full light-front calculation at NLO (full curve: scenario A, 
dotted curve: scenario B) is compared with the non-relativistic approximation (dot-dashed curve). 
CETQ4 fit of ref. [39]: full triangles. 
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FIG. 9. Left panel: unpolarized NLO (DIS) valence parton distribution at Q 2 = 10 GeV 2 as 
fuction of x. Light-front calculations: full curve; non-relativistic approximation (dot-dashed curve); 
NLO (DIS) CTEQ4 fit of ref. [39]: full triangles. Right panel: total unpolarized NLO (DIS) sea 
distribution at at Q 2 = 10 GeV 2 . Light-front results: scenario A (full curve), scenario B (dotted 
curve), non-relativistic approximation (dot-dashed curve). 
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